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Abstract The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) 
General Survey is a spectroscopic survey that will eventually cover approximately 
half of the celestial sphere and collect 10 million spectra of stars, galaxies and 
QSOs. Objects both in the pilot survey and the first year general survey are in¬ 
cluded in the LAMOST First Data Release (DR1). The pilot survey started in 
October 2011 and ended in June 2012, and the data have been released to the 
public as the LAMOST Pilot Data Release in August 2012. The general survey 
started in September 2012, and completed its first year of operation in June 2013. 
The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, 
of which 1,790,879 spectra have observed signal-to-noise ( S/N ) >10. All data with 
S/N >2 are formally released as LAMOST DR1 under the LAMOST data policy. 
This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stel¬ 
lar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not 
only spectra, but also three stellar catalogues with measured parameters: AFGK- 
type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 
entries), and M stars (121,522 entries). This paper introduces the survey design, 
the observational and instrumental limitations, data reduction and analysis, and 
some caveats. Description of the FITS structure of spectral files and parameter 
catalogues is also provided. 

Key words: techniques: spectroscopic survey — data release — catalog 


1 INTRODUCTION 

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) is a 4.0m telescope 
located at the Xinglong Observatory northeast of Beijing, China. The telescope is characterized 
by both a large field of view and large aperture. A total of 4000 fibers mounted on its focal plane 
make the LAMOST telescope a powerful tool with high spectral acquisition rate (Deng et ah, 
2012). With the capability of its large field of view and immense multiplexing, the LAMOST 
telescope is dedicated to a spectral survey of celestial objects over the entire available northern 
sky. After first light in 2008, commissioning observations were conducted for two years. In spite 
of the ongoing characterization of the LAMOST system, observation limitations (e.g., dome 
seeing, fiber positioning accuracy), and data processing software, the commissioning data have 
contributed some preliminary results on quasars (Huo et ah, 2010; Wu et ah, 2010a, b), planetary 
nebulae (Yuan et ah, 2010), and nearby galaxies (Zou et ah, 2011). 

In order to inspect the real instrumental performance and assess the feasibility of the science 
goals, the science working groups of LAMOST came up with the Pilot Survey, which was 
launched on 24 October 2011. In the Pilot Survey, both Galactic and extragalactic surveys were 
conducted - these were named LEGUE (LAMOST Experiment for Galactic Understanding 
and Exploration) and LEG AS (LAMOST Extra GAlactic Survey), respectively (Zhao et ah, 
2012). LEGUE is made up of three parts, which focus on the Galactic anti-center, the disk at 
longitudes away from the Galactic anti-center, and the halo. Some scientific results from the 
Pilot Survey include catalogs of white dwarf-main sequence binaries (Ren et ah, 2013), DA white 
dwarfs (Zhang et ah, 2013; Zhao et ah, 2013), and M dwarfs (Yi et ah, 2014), and identification 
of velocity substructure in the Galactic disk (Carlin et ah, 2013). Meanwhile, the LEG AS 
extragalactic portion of the survey studies galaxy redshifts, stellar populations of galaxies, 
and QSOs. Scientific results obtained include the identification of QSOs in the background of 
M31/M33 (Huo et ah, 2013), a study of the radio-loud quasar LAMOST J1131+3114 (Shi et ah, 
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2014a), and a search for double-peaked, narrow emission-line galaxies (Shi et ah, 2014b). The 
LAMOST Pilot Survey ended in June 2012, and included nine full moon cycles covering all 
good observing seasons at the site. The Pilot Survey data release includes 717,660 spectra, in 
which there are 648,820 stars, 2,723 galaxies, and 621 quasars, and a stellar parameter catalogue 
containing entries for 373,481 stars (Luo et ah, 2012). 

From the experience of the Pilot Survey, the LAMOST science committee fully analyzed the 
observations, assessing the magnitude limit, characteristics of the telescope, and data success 
rate. A five-year General Survey of LAMOST was designed, which would mainly focus on the 
Galactic survey of stars, yet still include a significant number of extra-galactic sources. 

As the only resolved classical disk galaxy, the Milky Way serves as a laboratory to study the 
formation and evolution of galaxies. However, there have been no large-scale Galactic spectro¬ 
scopic surveys with a uniform sampling over all stellar populations, due to the time consuming 
nature of observing a large fraction of its ^20 billion members spread over the entire sky. Thus 
our knowledge of the detailed structure of our Galaxy lags behind the understanding of the 
large-scale universe. Several notable achievements have been made by modern large-scale sur¬ 
veys such as Sloan Digital Sky Survey (SDSS); however, because SDSS is mainly focused on 
extragalactic studies, its spectroscopic observations have sampled only a small fraction of the 
Milky Way’s stellar components. The powerful spectral acquisition rate of LAMOST enables a 
homogeneous and (statistically) complete Galactic survey. 

The design of the survey will be shown in Section 2, including fields and sample selection. 
The observation conditions and limitations due the instruments will be discussed in Section 3. 
The LAMOST data reduction, spectral analysis, and stellar parameter pipelines will be intro¬ 
duced in Section 4. The data products of DR1 will be introduced in Section 5, as well as caveats. 
Finally, the data access and use are introduced in Section 6. 

2 DESIGN OF THE SURVEY 

2.1 LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) 

As a major component of the LAMOST project, the LEGUE spectroscopic survey of the 
Galactic halo will observe ^5.8 million objects with r < 16.8 mag at \b\ > 30° located mostly 
in the sky area with SDSS photometric coverage (Yang et ah, 2012; Deng et ah, 2012; Zhang 
et ah, 2012). Considering time consuming for faint objects which should be observed in dark 
nights, two medium-bright (M) plates stripes with high priority are proposed. One is the area 
at Dec. = 30° in Northern Galactic cap, and the other is at Dec. = 5° in Southern Galactic 
cap. 

Another unique component of LEGUE is the LAMOST Spectroscopic Survey of the Galactic 
Anti-center (LSS-GAC), which aims to survey a significant volume over a continuous sky area 
of ~ 3,400 deg 2 , covering Galactic longitudes 150° < t <210° and latitudes |6| < 30° (Yuan 
et ah, 2014a). The input catalog is selected from photometric catalogs of the Xuyi Schmidt 
Telescope Photometric Survey of the Galactic anti-center (XSTPS-GAC) (Liu et ah, 2014), 
with most targets selected with uniform probability. The catalog will include a statistically 
complete sample of about 3 million stars with magnitudes ranging from 14.0 < rl7.8 (and 
as deep as 18.5 mag for some fields), with rare objects of extreme colors targeted with high 
priorities. A detailed description of LSS-GAC is presented by Yuan et al. 2014a. 

For the Galactic disk survey, eight plates are chosen near the Galactic plane (|b| < 20°)and 
nearly uniformly distributed along the Galactic longitude, as well as in the region 0° < a < 67° 
and 42° < S < 59° (Chen et ah, 2012). This survey concentrates on open clusters (OCs), and 
is expected to obtain about 8.9 million stars down to a limiting magnitude of r ~ 16 mag. The 
highest observational priority is given to known or potential members of OCs, and a total of 
194 OCs have high priority. The remainder of the disk targets in the 11.3 < Imag < 16.3 range 
were uniformly selected over I and B—I color-magnitude space. Very high proper motion objects 
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(7mas/yr) are removed from the target catalogs. Positions, proper motions, and magnitudes are 
taken from the PPMXL (Roeser et ah, 2010) catalog. 

2.2 LEGAS (LAMOST ExtraGAlactic Survey) 

The extra-galactic survey (LEGAS) consists of two main parts - a galaxy survey and a QSO 
survey. For galaxies, it is split into two regions (again, to avoid the Galactic plane). One is 
in the North Galactic Cap region, and shares the same footprint as the SDSS legacy ellipse. 
This survey mainly aims to observe those SDSS main sample galaxies (with r < 17.75) that 
were not spectroscopically observed by SDSS due to fiber collisions (Stoughton et ah, 2002). 
Such “missed” SDSS galaxies are named as the complementary galaxy sample in the LAMOST 
survey. There are 68,722 such targets, which are given highest priority for fiber assignment 
on plates when possible. Because the number density of galaxies is very low, the remaining 
fibers on LEGAS plates are filled by LEGUE halo survey stars. The other region of the galaxy 
survey in LAMOST is the South Galactic Cap region (SGC). The LEGAS team has a long-term 
strategic goal of a galaxy survey in the SGC footprint of LAMOST (b < —30° and S > —10°), 
which aims to take spectra of all galaxies with r < 18, and a sample of blue galaxies down to 
r < 18.8. Since this desired galaxy survey in the SGC requires much deeper observations than 
the main LAMOST stellar sample (currently, mostly limited to r < 16.8), only a few dark nights 
per observing season are available for these galaxies. Therefore, the galaxy survey in SGC now 
focuses on a strip, ranging from 45° < a < 60° and 0.5° < S < 9.5°. Besides these two well- 
defined samples of galaxies, some bright-infrared galaxies are selected as extra observational 
targets which are selected from infrared surveys such as IRAS, WISE, and HERSCHEL. 

The efficiency of identifying quasars with redshifts between 2.2 < 2 < 3 is low in SDSS 
(Schneider et ah, 2010), because quasars with such redshifts usually have similar optical colors 
as stars and are thus mostly ignored by the SDSS quasar candidate selection algorithm. As 
permitted by the observational conditions and instrument standards of the LAMOST telescope, 
LAMOST quasar survey aims to find as many quasars as possible. We selected quasar targets 
from optical data within the Sloan photometric footprint using data mining methods such as 
support vector machines (SVM) (Peng et ah, 2012; Zhang & Zhao, 2003), distribution of quasars 
and stars in different color spaces including Y — K vs. g — z (Wu et ah, 2012)) and z — W1 
vs. g — z (Wu & Jia, 2010), and the extreme-deconvolution method (Bovy et ah, 2011). The 
targets selected by multiple methods are given much higher priority in observation. The core 
sample is contributed from the above methods using survey data from SDSS, UKIDSS, WISE. 
The bonus sample comes from other survey data, such as X-ray bands from ROSAT, XMM, 
and Chandra, and radio bands from FIRST and NVSS. It is noted that the QSO targets are 
not limited to the LEGAS plates, and are given high priority in the LEGUE plates as well. 

2.3 The special fields 

2.3.1 M31-M3 3 

Targets in the M31-M33 region include known and candidate planetary nebulae (PNe), HII 
regions, super-giants, and globular clusters in M31 and M33, plus known and candidate back¬ 
ground quasars and foreground Galactic stars selected from the XSTPS-GAC photometric cat¬ 
alogs. These stars are selected, assigned and observed in the same way as for the LSS-GAC 
main survey. The selection criteria of various types of targets and candidates of interest are 
presented elsewhere (e.g., Yuan et al. 2010 for PNe; Huo et al. 2010 for quasars). The field 
central stars for the M31-M33 survey are selected similarly to the LSS-GAC main survey in 
the (RA, Dec) plane. Two groups of field central stars are chosen to cover the M31-M33 area 
(0° < RA < 30°, 25° < Dec < 50°) independently. Each group contains 41 central stars. 
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2.3.2 The LAMOST Complete Spectroscopic Survey of Pointing Area at Southern Galactic 
Cap 

A LAMOST Key Project named the LAMOST Complete Spectroscopic Survey of Pointing 
Area (LCSSPA) at SGC is being conducted (Yang et ah, 2015), which is designed to have 
repeated spectroscopic observations of all sources (Galactic and extra-galactic) in two 20 deg 2 
fields in the SGC. The central coordinates of the fields are (a, S) = (37.88150939°, 3.43934500°) 
and (21.525988792°, —2.200949833°), respectively. Targets in these field mainly consist of stars, 
galaxies, ix-band variables, and QSOs. The stars and galaxies are selected from SDSS Data 
Release Nine (DR9) using r-band psf magnitudes and petro mags between 14.0 < r < 18.1 
mag, respectively. The riband variable sources are selected as objects for which the difference 
in magnitudes between the SDSS and South Galactic Cap r^-band Sky Survey (SCUSS) are 
larger than 0.2 mag; these are a magnitude-limited sample down to u = 19.0 mag. The QSOs 
are selected by the same methods discussed above. The highest priority is assigned to riband 
variables and QSOs, with galaxies are the second priority level, and stars having the lowest 
priority. 

2.3.3 Kepler fields 

Additionally, a total of 14 Kepler fields with a fixed field of view (FoV) of 105 square degrees 
in the constellation Lyra and Cygnus are also selected according to the prioritized list provided 
by LK-project (Cat et ah, 2015). The prioritized lists are designed mainly based on the type 
of targets and brightness. According to the object types, these targets are divided into four 
priority levels: standard targets (~250, MK secondary standard stars), KASC targets (^7,000, 
selected by the Kepler Asteroseismic Science Consortium, or KASC), planet targets (^150,000, 
selected by the Kepler planet search group Batalha et al. 2010), and extra targets (^1,000,000 
other targets in the Kepler FoV from the KIC 10). For these, the higher observation priority is 
given to the fainter objects. 

2.4 Plate design 

2.4.1 Plate definition and their magnitude 

The LAMOST plates are divided into four modes, very bright (VB), bright (B), medium-bright 
(M), and faint (F) according the magnitude range of stars to be observed. In DR1, the VB mode 
contains 2693 plates, and is used to observe stars brighter than r = 14 mag over the entire sky 
observable by LAMOST (-10° < Dec <60° ) in order to make use of bright nights or nights of 
unfavorable conditions (e.g., large seeing or low atmospheric transparency). In the XSTPS-GAC 
footprint, stars are selected as potential targets from the XSTPS-GAC and 2MASS catalogs 
with magnitudes r < 14.0 mag and 9.0 < J < 12.5, respectively. Stars outside of the XSTPS- 
GAC sky area are selected as potential targets from the PPMXL (Roeser et ah, 2010) catalog 
with magnitude 10.0 < bl < 15.0, 10.0 < b2 < 15.0, 9.0 < rl < 14.0, 9.0 < r2 < 14.0, 
or 8.5 < i < 13.5, and from the 2MASS catalog with magnitude 9.0 < J < 12.5. The 
“B” (bright) observing mode contains 751 plates in DR1 that observed stars with r magnitudes 
between between 14 and 16.8. There are 692 and 655 plates assigned to the M (16.8 < r < 17.8), 
and F (r > 17.8) designations, respectively. 

2.4-2 Tiling method (Cui et al, 2012) 

Merging all of the science goals described above, the science committee of LAMOST designed 
the survey strategies and input catalogues. Based on these, the LAMOST team developed a 
software named SSS (Sky Strategy System) to produce observational plates that most efficiently 
use the multi-fiber abilities of LAMOST. The maximum density algorithm, density gradient 
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algorithm, and the mean-shift algorithm are used for most sky areas, with the exception of the 
Galactic Anti-Center Survey, which uses separate target selection criteria (Yuan et ah, 2014a). 

Given a sky area, the maximum density method can find the tile (or plate) with the maxi¬ 
mum object count by making a uniform tiling. By traversing all the tiles covering the sky area, 
a tile with the maximum object count will be found and then a new observational plate will 
be built at that position. The only factor that affects the speed and accuracy of this algorithm 
is the target distribution in the sky area. The denser the covering is, the higher the accuracy 
of finding the maximum tile becomes and the more efficient is the method, and vice-versa. 
The density gradient method is essentially iterative. Based on a uniform covering, this method 
makes an estimate of the object distribution and moves the tile in the direction that is expected 
to have a higher object density. The basic idea of the mean-shift algorithm is also iterative. As 
long as the objects are distributed unevenly in the sky areas, a density gradient exists between 
different positions. Vectors are used to describe the position of targets in the tile, and the sum 
of vectors is set along the direction of the density gradient. By iteratively moving the field 
center along the sum of vectors and calculating the new vector summation, the recent peak 
density location is obtained. In order to maintain the goal of uniformity of the remaining area 
and to avoid neglecting some of the peak positions in areas where density varies intricately, the 
algorithm will try to find a position with a global higher density and make a quick coverage on 
the whole available sky area before iteration (Ren et ah, 2007). 

For each plate, one bright star (the so called central star, brighter than 8 mag) must be 
placed in the plate center to feed a Shack-Hartmann detector that drives the active optics 
system, and four stars (brighter than 17 mag) must be placed on each of four guiding CCD 
cameras. The bright stars located in the center were selected from Hipparcos catalog (Perryman 
et ah, 1997), with binaries excluded, and at the time of observation, the focal plane is rotated a 
bit to confirm that the proper guiding stars are inside the FoV of the four guiding CCD cameras. 
Guide stars were selected from the GSC2.3 catalog (Lasker et ah, 2008). Central stars will be 
excluded if there are no corresponding guide stars. Since in some areas (e.g., at low Galactic 
latitudes) there are many more available central stars than needed to fully cover that part of 
sky, a step by step search was applied to reject redundant bright stars and find the tiling that 
maximizes sky coverage of that area (Yuan et ah, 2012). 

2.4-3 Fiber allocation scheme 

Once the central position of a plate is decided by the tilling method, each object in the input 
catalog will be assigned a priority according to its scientific requirements. Five flux standards 
and 20 sky fibers are reserved in each spectrograph for data reduction; therefore, these five flux 
standards are allocated with the highest priority. All of the scientific targets are assigned in 
order of priority; however, if more than two targets have the same priority in the same fiber 
cell, the closest target will be chosen to avoid collision with other fibers. On the other hand, 
if a target can be reached by different fibers at the same time, then it will be assigned to the 
fiber within its fiber cell that contains the minimum number of possible targets (Yuan et ah, 
2014b). A total of 20 sky fibers will eventually be assigned from fibers for which no targets 
were assigned. If the number of sky fibers allocated to a given spectrograph is less than 20, the 
lowest priority fibers will be reassigned to the sky until there are 20 sky fibers. 

3 LIMITATIONS ON THE OBSERVATIONS AND PERFORMANCE OF THE 
TELESCOPE 

3.1 Seeing 

One of the significant factors affecting the observations is seeing, which consists of the native 
site seeing as well as dome seeing. The seeing depends on the season, and is worse in the windy 
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winter, but with a huge diversity. The peak in the distribution of the nightly site seeing is 
at around 3 arcsec, and 85% of the seeing measurements are better than 4 arcsec (Yao et ah, 
2012). The non-classical dome of the LAMOST telescope seemingly affects the seeing measured 
directly at the focal plane, with the measurements showing seeing of about 3.5 arcsec on average. 
Compared to the site seeing, it is dome seeing that contributes the most in limitating the system 
capabilities (especially affecting the limiting magnitude). Considering the fiber diameter of 
3 arcsec, an average increase of 1 arcsec in seeing (at the focal plane) will lead to a 0.4-magnitude 
decrease in flux collected by each fiber. 


3.2 Fiber positioning and use efficiency 

The precision of fiber positioning is an important factor affecting the efficiency of the survey. 
Through repeated observations and calculations, the accuracy of fiber positioning has been 
determined to be less than 1.5 arcsec on average (Chen et ah, 2014), which has less effect on 
the fiber magnitude than the seeing. Another limitation is provided by the fixed fiber density 
of around 200 per deg 2 . This limits sampling in both dense and sparse fields in different sky 
areas. For dense fields, objects are too crowded to resolve (especially near objects such as 
stellar clusters). In contrast, the fiber density results in a considerable waste of fibers (which 
are allocated to sky) for sparse fields in high Galactic latitude areas. 


3.3 The performance of the instruments 

The performance of the instruments plays the same important role in the quantity and quality 
of spectra observed as the procedure of data reductions. We briefly describe it as follows. 


Overhead During the observations, the overhead time will reduce the observational efficiency 
of the survey to some extent. Specifically, the telescope needs to be confocal, guide stars need to 
be located, and the fibers need to be repositioned when a new plate is to be observed. Including 
the CCD readout time, the total overhead prior to observation of a field takes about 20 minutes. 

Throughput The throughput is a major indicator we use to measure the efficiency of 
the telescope. The whole throughput takes the contributions of the spectrograph, fibers, CCDs, 
optical system, etc., into account. In general, the throughput of the red and the blue channels 
are 4% and 2% respectively Cui et al. (2012). During the pilot survey and the 1st year regular 
survey, the throughput of LAMOST was relatively low, and also showed large dispersion among 
the 4,000 fibers. The 4 th spectrograph (out of a total of 16 spectrographs) has the highest 
throughput of those located near the center of focal plane; the median throughput of the valid 
fibers connected to the 4f h spectrograph during the pilot survey was about 1.0% and 2.0% for 
the blue and red bands, and improved to 1.2% and 3.0% for the blue and red bands during the 
1st year regular survey. The No. 3, 5, 8, 9, and 15 spectrographs are located at the middle-ring 
of the focal plane; their throughputs are about 80—90% of that of the 4 th spectrograph. The 
No. 1, 2, 6, 7, 10, 11, 12, 13, 14, and 16 spectrographs are located at the outer-ring of the focal 
plane, and have measured throughputs that are about 50—80% of that of the A th spectrograph. 
Among the 250 fibers accommodated by each spectrograph, the difference in the throughput 
can change by a factor of two. 

Stability of instrument profile Variations of emission lines from arc lamps have been 
monitored during a continuous 300 day period in order to measure the stability of the 
instrument profile. From Figure 1, it can be seen that the profile of each fiber changes little 
with time, but that there differences among the profiles between the fibers. 
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Fig. 1 Variations of arc lines during continuous 300 days. The left panels show Full- 
Width-Half-Maximum (FWHM) of arc lines in the blue channel for the first 100, 
the second 100 and the last 100 days respectively, and the right panels are the same 
observations for the red channel. 


Flux calibration. The LAMOST project is designed as a spectroscopic survey without a 
photometry telescope equipped. Therefore, absolute flux calibration cannot be done without the 
photometric data (Song et ah, 2012). Besides, not enough first-level standard stars are available 
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over the large FoV observed by LAMOST that can be used for relative flux calibration. Thus 
we choose some stars with high quality spectra as standard stars. The uncertainty of the shape 
of the continuum in these standard stars will thus contribute uncertainty to the LAMOST flux 
calibration. 

3.4 Sky subtraction and flat fielding 

According to the requirements set for LAMOST observations, ^ 10 — 15% of the fibers for each 
plate are assigned for modelling the sky. This fiber assignment strategy works well for fields 
at high Galactic latitude, providing sky-subtraction that is consistent with results from the 
SDSS survey. However, the proportion of the fibers that sample “clean” sky background is not 
adequate to build a supersky in observations near the Galactic disk, especially in bright nights. 
The emission lines of the night sky spectra are not subtracted satisfactorily for a number 
of spectra under these conditions. Meanwhile, the fiber-to-fiber profile differences makes the 
subtraction of sky emission lines more complicated. Moreover, the sky continua are more difficult 
to remove than narrow sky emission lines, such that there will always be at least some sky 
residual remaining in reduced spectra. This in turn presents a challenge we have to deal with 
when determining (stellar) parameters from spectra. 

It is not easy to find a light source with a five-degree FoV to use as a flat field. We therefore 
separately correct the flat field for each spectrograph, since each spectrograph covers only about 
a 1° FoV. Twilight flats and screen flats are used. Additionally, the vignetting and the change 
of aperture caused by the telescope pointing also complicate the correction of the flat field. The 
movement of fibers may lead to focal ratio degradation caused by fiber tension changes, which 
will cause differences in the efficiency of fibers. 

4 DATA REDUCTION 
4.1 Statistics of data quality 

4-1.1 Resolution and resampling 

The theoretical resolution of the 16 LAMOST spectrographs is R = 1000 (as determined by 
the gratings), and the designed wavelength coverage is 365-900 nm. In practice, the resolution 
reaches R>1500 by narrowing the slit at the output end of the fibers to 2/3 the fiber width. Since 
the spectral images are sampled on CCDs, to avoid confusion between resolution and sampling 
for the final re-binned spectrum, we explain the difference as follows. On the original CCD 
image, 4K pixels are used for recording blue or red wavelength regions ranging from 370 nm to 
590 nm or from 570 nm to 900 nm respectively, which means each A has been sampled in 2 CCD 
pixels in the raw data. When the blue and red channels are combined together, each spectrum 
is re-binned for easier calculation of radial velocity. The combined spectra are re-sampled in 
constant-velocity pixels, with a pixel scale of 69 km s -1 , which means the wavelength difference 
between two adjacent points is A log (A) = 0.0001. For convenience, the re-binned data interval 
is still referred to as ’pixel’. 

4.1.2 The definition of Signal to Noise Ratio 

The signal to noise ratio (SNR) is always used as a quality indicator, and is the ratio of flux to 
noise in some bandpass. Generally, SNR is an average value in a wavelength band and represents 
predominantly the continuum flux level. SNR hereafter indicates the SNR per ’pixel’. In the 
LAMOST data products, SNR is defined by using the concept of inverse-variance, which is 
included in the spectral FITS files. It is easy to obtain the SNR for each pixel by calculating 
the product between the flux and its inverse-variance. The mean SNR in a wavelength band can 
be averaged from the SNR of each pixel in the band range. Another way to obtain SNR is to 
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Fig. 2 The two definitions of SNR are compared for a sample of L AMO ST spectra. 
From left to right, the panels show SNR in the g , r, and i bands. 


Table 1 The statistics of spectrum success rate 


Plate type 

Number 

of plates 

Number of 

spectrographs 

Rate of 

successful spectra 

GAC 

190 

2592 

318854/531019 

= 60% 

M31 

41 

511 

54950/104770 

= 52% 

VB 

343 

4405 

552245/794100 

= 70% 

HD 

347 

5124 

479076/745791 

= 64% 

EG 

71 

1057 

102675/199091 

= 52% 

F M B 

202 

2805 

268154/570498 

= 47% 

Kepler 

8 

104 

15333/22899 

= 67% 

Total 

1202 

16598 

1791287/2968168 

= 60% 


Notes: 

1. The last column is the spectrum success rate, which means the ratio of the number of spectra with 
SNR^ or SNR r > 10 to the number of fibers that were assigned to targets. 

2. Plate types are named from the pilot survey and the first year of the general survey. The details are 
described as follows: 

GAC : Galactic Anti-Center plates. 

VB : Very bright stars with r<14 in all northern sky. 

HD: Halo and disk plates (Sep. 2012—Jun. 2013). 

EG: Extragalactic plates. 

F M B: Halo and disk plates (Oct. 2011—Jun. 2012). F, M, and B means faint, medium, and 
bright, respectively. 


calculate the ratio between the continuum and the la standard variance. These two definitions 
of SNR are shown to provide consistent results when applied to LAMOST data (Figure 2). 

4-1.3 Data quality 

The data quality is determined by SNR, and a successful observational target should have SNR 
>10 in g— or i— band. Table 4.1.3 lists the statistics of the spectrum success rate. The success 
rate improved from the pilot survey to the general survey, and the rate on VB plates reaches 
to 70% even though these targets were observed in bright nights. 

4.2 2D pipeline 

Raw CCD data were reduced by the LAMOST data reduction software known as the LAMOST 
2D pipeline. The procedures of the 2D pipeline, similar to those of SDSS (Stoughton et ah, 
2002), aim to extract spectra from CCD images and calibrate them. The main tasks of the 
2D pipeline include dark and bias subtraction, flat field correction, spectra extraction, sky 
subtraction, wavelength calibration, sub-exposure merging and wavelength band combination. 
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The processes of the 2D pipeline are listed as following (see the flowchart in Figure 3). 

1. Bias and dark counts are subtracted from each raw image. During the pilot survey, some 
inside CCD light source affected images, so that we needed the dark frames to eliminate this 
stray light. In the general survey, this light source was identified and removed, so that no dark 
subtraction is needed. 

2. Flat field corrections are mainly achieved using twilight flat fields. However, white-screen 
flat fields were used on some days. LAMOST spectrographs are regarded as 16 different 
individual instrument systems for flat correction, so that the flat for each spectrograph is 
independent. 

3. The flat-field spectra are traced for each fiber, and the centroid of the row position of 
each fiber is fitted by a polynomial function. A good fitting of each spectral track overcomes 
interference between fibers. However, the fitting will fail at the end of the wavelength coverage 
range when the SNR becomes low. 

4. A spectrum is extracted along the track fitted by step 3 using a weighted aperture 
integration method; flux on 15 pixels is integrated to provide a single flux for ID spectra. 

5. The sky background is modeled and subtracted by a fixed proportion of the fibers on 
each plate that is reserved for pointing at blank sky. The method of B-spline fitting is applied 
to make a supersky for each sub-exposure and then Principal Component Analysis (PCA) is 
used to recorrect the supersky after combining in this procedure. Additionally, the relative 
strength of sky emission lines is used to correct efficiency differences between sky fibers and 
object fibers, since these lines exists both in sky spectra and object spectra. 

6. Wavelength calibration is by use of the arc lines. The arc lamp spectra are extracted, and 
centroids of the lines are measured, to which we fit a Legendre polynomial as a function between 
wavelengths and pixels, with a fifth-order polynomial in the blue branch and sixth-order in the 
red. A vacuum wavelength scale is applied for wavelength calibration, is slightly adjusted to 
match the known positions of certain sky lines, and then corrected to the heliocentric frame. 

7. The telluric absorption in four wavelength regions (6850—6960 A, 7150—7350 A, 
7560—7720 A, 8100—8240 A) is removed. 

8. The spectra are flux-calibrated by matching the selected flux standard stars with 
their matching templates of F-type stars. First, 250 spectra from one spectrograph are 
sorted in descending order of SNR. Stars with bad SNR or saturated fibers are abandoned. 
After normalizing the continuum of each spectrum, atmospheric parameters are obtained by 
matching a template using the Kurucz model. Then stars with 6000 K < T e ^ < 7000 K are 
selected as standard stars, after visual inspection to confirm their correct identification. If the 
number of standard stars is less than 3, then the effective temperature will be extended from 
5500-7500 K to obtain more standard stars. 

9. For each object in all the exposures, the spectra from the red and blue spectrograph are 
combined by stacking the points with corresponding wavelengths using a B-spline function, 
with inverse-variance weighting. Outliers due to cosmic rays are removed and masked, then 
errors of the fluxes are estimated. 
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Fig. 3 Flowchart illustrating the LAMOST 2D pipeline processing. 


4-2.1 Uncertainties in the wavelength calibration 

During the wavelength calibration, errors are most likely to be introduced when finding centroids 
of the emission lines of the arc lamp. In these emission lines, not all of the lines are single - 
some observed lines are actually blends of two or three lines. Therefore, the fitting quality of 
a unimodal function for these mixtures may have large uncertainties. Statistically, the average 
error of wavelength calibration for the arc lamp is less than 0.02 A , corresponding to about 
1 km s _1 . Generally, the calibration accuracy of blue spectra is better than that of the red one, 
but it degrades below 4000 A. 

Another error of wavelength calibration is from heliocentric correction by wrongly using 
Beijing local MJM instead of International Standard Time. The difference between local MJM 
and International Standard Time is 480 minutes, which results in a random error of 0.36 km 


4-2.2 Uncertainties in the relative flux calibration 

If there are no standard stars observed on one of the spectrographs (for a given plate), combined 
object spectra spanning the entire wavelength range cannot be produced. This is because the 
matching of the blue- and red-channel spectra relies on matching their flux in the wavelength 
region where they overlap. 

The reddening calibration is not considered in standard stars selection. The extinction is 
negligible in high Galactic latitude areas. However, reddening will affect flux calibration severely 
in the low Galactic latitude areas, biasing the continuum shapes of resulting spectra in these 
high-extinction regions. 

The flux of calibrated spectra from different exposures for each object may be not on the 
same scale. Hence before combining them, we adjust them onto a uniform scale, and then 
combine the spectra from different exposures by the method of B-spline fitting, which is similar 
to the SDSS pipeline. 

4.3 ID pipeline 

The LAMOST ID pipeline works on spectral type classification and measurement of radial 
velocity (RV) for stars or redshift for galaxies and QSOs. By using a cross-correlation method, 
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the pipeline recognizes the spectral classes and simultaneously determines the redshifts or radial 
velocities from the best fit correlation function. In this procedure, a low order polynomial is 
added while matching the observed spectra with templates, since the physical continuum can 
be difficult to detect because of the low resolution, reddening, and uncertainties in the flux 
calibration. Currently we employ the same templates as SDSS. However, the templates will be 
updated with LAMOST spectra for reductions after DR1. Finally, the pipeline produces four 
primary classifications, namely STAR, GALAXY, QSO, and UNKNOWN. 

Measuring redshifts and recognizing galaxy and QSO spectra are difficult tasks since the 
pipeline does not work as well as for stellar classification. An additional independent pipeline 
was designed for galaxy recognition and redshift measurement after the initial ID pipeline is 
run. It automatically identifies galaxies from the objects with confidence of classification below 
80 or a target type as “galaxy” in the input catalogue. The redshifts of galaxies are measured 
through line centers; a Gaussian function with sigma of 1.5 times the wavelength step is applied 
to the spectra to eliminate noise before line centers are measured. The continua smoothed by 
a median filter are divided out to complete the normalization. Those data points that exceed 
2cr of a normalized spectrum are selected as candidate emission lines, and a set of Gaussian 
functions are used to fit the lines. All the line centers are matched with line lists spaced by steps 
of 0.0005 in redshift (z). If most of the lines are matched successfully with heavily weighted 
lines such as H a , Oil, Hp, OIII, or Nil for emission line galaxies, or NaD, Mgb, Call H, or Call 
K for absorption galaxies, the spectrum is classified as galaxy, and the corresponding z is the 
raw redshift of the spectrum. To detect mistakes in the redshift measurement, this redshift is 
compared with the results of the ID pipeline, and the final redshifts are decided by experts. For 
QSOs, the recognitions and measurements highly depend on visual inspection; a QSO catalogue 
of DR1 detections will be released by Ai et al. (2015). 

4.4 LAMOST Stellar Parameter Pipeline - LASP 

LAMOST has carried out the pilot survey and the first year general survey, producing more 
than 2.2 million spectra, among which nearly 2 million were stellar spectra. Through statistical 
analysis of the stellar physical properties of this tremendous spectroscopic database spanning 
wide sky coverage, LAMOST will definitely encourage and actively promote our understanding 
on the structure, chemistry, and kinematics of the Galactic disk and the halo as well as for 
its origin and evolution history. Efficient parameterization software to accurately derive the 
fundamental stellar atmospheric parameters plays a key role for the scientific exploration. For 
this purpose, we developed and implemented the LAMOST stellar parameter pipeline (LASP, 
version 1.0). 

4-4-1 Input Data and Strategy for LASP 

LASP has been used to automatically determine the stellar parameters (effective temperature: 
T e ff, surface gravity: logg, metallicity: [Fe/H], and radial velocity: V r ) for the AFGK-type 
stellar observations from the LAMOST survey. The input data for LASP are the output FITS 
spectra produced by the 2D and ID pipeline. For the DR1 data release, the fits file data version 
is v2.6.4 (version according to 2D and ID pipeline updating status, for DR2 data release the 
version will be upgraded to v2.7.5 after some improvement), with a resolving power of R ~ 1800 
covering wavelengths 3800—9000 A. The survey observing nights are generally divided into three 
categories: dark nights (8 days before and after the new moon), test nights (3 days around the 
full moon, used for instrumental testing), and bright nights (the remaining nights in a lunar 
cycle). The current criteria required for stars to be processed by the LASP are: final .class is 
STAR, and final.subclass is spectral type A, F, G, or K, and g-band SNR of SNR g > 15 and 
SNRg > 6 for the bright and dark nights, respectively. The final.class and final.subclass 
information is retrieved from the ID pipeline results. 
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4-4-2 Methods of the LASP 

The LASP consecutively adopts two methods - CFI and UlySS Wu et al. 2011 b - to determine 
the stellar parameters. The CFI method is used to give a set of coarse measurements that serve 
as initial guesses for ULySS, which then gives the final measurement that is issued by LASP. 

1. Correlation Function Interpolation (CFI) method: Initial stellar parameter estimates are 

provided by a method called the CFI, short for Correlation Function Interpolation (Du et al., 
2012 ). The core algorithm of the CFI method is as follows: provided the observed flux vector 
is O, and the synthetic flux vector is S', the best-fit case is cos < O, S > = 1 theoretically. 
The algorithm searches for the best-fit by maximizing cos < 0,S > as functions of T e ff,logg, 
and [Fe/H], where cos < O, S > = which is also referred to as correlation coefficient. 

Practically, CFI first explores in temperature among the grid of synthetic spectra to determine 
the T e ff value of the input spectrum. After T e # is fixed, the code searches for the [Fe/H] in 
metallicity space among the grid, and then with fixed T e ff and [Fe/H], it searches for the \ogg. 
According to maps of the correlation coefficient distribution for the 3 parameters, the reliability 
for T e ff estimation is relatively more credible, the [Fe/H] is less credible, and it is hard for the 
algorithm to determine a log# solution. Thus the biweight mean of several log# estimations 
with top correlation coefficient values is given out as the output. 

For the production of the synthetic spectra grid used by our CFI algorithm, we em¬ 
ployed Kurucz spectrum synthesis code based on the ATLAS9 stellar atmosphere models 
provided by Castelli & Kurucz (2003). This grid is smoothed to the resolution of the 
LAMOST spectrograph (R ~ 1800) and the pseudo-continuum of the grid is manually 
obtained by an experienced spectral analyser. The analyser manually marks the points 
which may be on the pseudo-continuum based on spectral features. Then spline curves are 
employed to fit the spectral pseudo-continuum. Finally, the points can be adjusted until 
the pseudo-continuum meet the requirements. The CFI’s adopted synthetic library contains 
8903 spectra, in a grid with stellar parameter coverage of: 3600iF < T e ^ < 7500 K in steps of 
200 K, O.Odex < log <7 < 5.0 dex in steps of 0.25 dex, and —3.0 < [Fe/H] < 0.4 in steps of 0.2 dex. 

2. ULySS method: The ULySS (Universite de Lyon Spectroscopic analysis Software, (Koleva 
et al., 2009; Wu et al., 2011b), available at http : //ulyss .univ-lyonl. fr/ package is employed 
by LASP to derive the stellar parameters via minimizing the squared difference between the 
observations and the model. The model is adjusted at the same resolution and sampling as the 
observation, and the fit is performed in the pixel space. This method determines all free param¬ 
eters within a single fit in order to properly handle the degeneracy between some parameters, 
e.g., temperature and metallicity. Wu et al. (2011b) have already described the y 2 minimization 
of the ULySS. The parametric model is: 


Obs(X) = P n (X) x [TGM(T ef f, log#, [Fe/H], A) ® G(v sys: a)] 

where A is the logarithm of the wavelength, Obs(X) is the observed stellar spectrum sampled 
in log A, P n (A) a series of Legendre polynomials of degree n, and G(v sys ,cr) is a Gaussian broad¬ 
ening function characterized by residual velocity v sys , and the dispersion a. The a encompasses 
both the instrumental broadening and the effects of rotation.The multiplicative polynomial is 
used to absorb errors in the flux calibration, Galactic extinction or any other source affecting 
the shape of the spectrum. It replaces the prior rectification or normalization to the pseudo¬ 
continuum that other methods require. The parameters of the TGM are evaluated together 
with those of the LOSVD with a Levenberg-Marquardt (Marquardt, 1963; More et al., 1980) 
routine (hereafter LM). The (linear) coefficients of the P n (A) polynomials are determined by 
ordinary least-squares at each evaluation of the function minimized by the LM routine. The 
gaps and bad pixels (e.g. emission lines, cosmic ray, bad sky subtraction, telluric lines or bad 
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ones due to the instrument errors etc.) are automatically rejected during the fitting process by 
ULySS, using the clipping algorithm applied iteratively on the residuals to the fit. 

The non-linear TGM function is an interpolator that can return a spectrum for a given 
set of atmospheric parameters [T e ff,log#, [Fe/H]] by making an interpolation over the whole 
ELODIE library (Prugniel & Soubiran, 2001, 2004; Prugniel et ah, 2007). The interpolator 
(available at http://ulyss.univ-lyonl.fr/models.html) consists of polynomial expansions 
of each wavelength element in powers of logT e ff, log#, [Fe/H], and f(a) (a function of the 
rotational broadening parameterized by cr, the standard deviation of a Gaussian). Three sets 
of polynomials are defined for three temperature ranges (roughly matching OBA, FGK, and M 
types) with important overlap between each other where they are linearly interpolated. For the 
FGK and M polynomials, 26 terms are adopted; for OBA, 19 terms are used. The coefficients 
of these polynomials were fitted over the ~ 2000 spectra of the library, and the choice of 
the terms of the T e ff limits and of weights were fine tuned to minimize the residuals between 
the observations and the interpolated spectra. This interpolator built on the ELODIE library 
(version 3.2) provides valid inverted atmospheric parameters covering 3100 ^ 59000 K in T e ff, 
0.00 ~ 5.00 dex in log#, and —2.80 ~ 1.00 dex in [Fe/H], from O to M type stars. The derived 
intrinsic external accuracies of ULySS for high quality AFGK stellar spectra are 43 K, 0.13 dex, 
and 0.05 dex for T e ff, log#, and [Fe/H] respectively (Wu et ah, 2011b). The performances of 
ULySS applied to LAMOST commissioning period observations were already illustrated in Wu 
et al. (2011a). 

4-4-3 Processing Procedure of the LASP 

LASP operates in two stages to measure the stellar parameters. Since the LAMOST spectral 
flux calibration is relative (because there is no valid complete system for photometric standard 
stars), in the first stage we measure the original spectrum, and in the second stage we use 
the normalized spectrum. For the medium resolution of the LAMOST survey spectra, the blue 
band spectra contain most of the physical features necessary to constrain the stellar parameters. 
To avoid the effects of low instrumental response near both wavelength edges, as well as for 
optimizing the computing time and using minimal storage space, CFI selects wavelengths of 
3850 — 5500 A, and ULySS selects 4400 — 5700 A(in the 1st stage) and 4400 — 6800 A(2nd stage) 
for the fitting range. 

In each stage, we first utilize CFI to get a set of initial coarse estimates. After this, we use 
CFI results as starting guesses for ULySS to obtain more accurate and credible measurements. 
Figure 4 displays a general flowchart of the LASP. 

As CFI runs its core algorithm to determine the parameters individually one after another, 
it determines the best match between the observation and the synthetic grid spectra, both of 
which have been normalized. In the first stage of CFI, it uses the radial velocity (RV) estimated 
by the ID pipeline to shift the observed spectrum to the rest frame. Theoretically, observed 
spectra could be approximated by synthetic spectra multiplied by low order polynomials. To 
quickly get an initial guess of the effective temperature for CFI as a prior, the code uses the 
NGS1 (Lee et ah, 2008) grid. The NGS1 grid of synthetic spectra was generated from the 
Kurucz NEWODF models, which employ solar relative abundances from Grevesse & Sauval 
(1998). More details about generation of NGS1 synthetic spectra were described in Lee et al. 
(2008). These synthetic spectra were convolved to R = 1000, then sampled into 1.67 A per 
pixel to be the NGS1 grid spectra. In order to speed up the analysis, only the spectral range of 
4500 A to 5500 A was used. The synthetic spectra have been normalized by using a 9th order 
polynomial, iteratively rejecting points that are more than la below and 4cr above the fitted 
function. The NGS1 grid covers 4000 K<T e ff <9750 K in steps of 250 K, 0.0 dex< log g <5.0 dex 
in steps of 0.25 dex, and —4.0 < [Fe/H] < +0.5 dex in steps of 0.25 dex for the three parameters 
respectively. The observed LAMOST spectra are also processed in the same fashion over the 
same wavelength range, then we search for the best-fit using the CFI core algorithm to get an 
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Fig. 4 Flow chart illustrating the parameter estimation process of LASP. 


initial guess of T e ^. Then with this T e ff, CFI uses its core algorithm to search for the most 
similar synthetic spectra (with already known pseudo-continuum, which could be subtracted 
before fitting) in the grid described before in Sec.4.4.2. We use a 5th order polynomial to absorb 
the difference between the normalized observation and the normalized synthetic spectra, and 
finally get an estimation of the three parameters one after another. The grid described in Sec. 
4.4.2 is more intensive than the one of NGS1. In the second stage of CFI, it takes the ULySS 
first stage results [T e ff, logg, [Fe/H], RV], which provide a more accurate RV to correct to the 
rest frame, and removes the need to use the NGS1 grid, since we can use the T e ff given by the 
ULySS 1st result. Since the spectra dealt with in the second stage are normalized, CFI directly 
matches the spectra without a prior step removing the pseudo-continuum. Figure 5 shows an 
example of the CFI spectral fitting both in the blue range and whole range. 

Before stellar spectra are parameterized, the pseudo-continuum must be divided out (or 
normalized). For the input observation, its analogous synthetic spectrum can be found in the 
synthetic grid (Sec 4.4.2) according to the parameters of ULySS 1st stage results. Since the 
pseudo-continuum of the grid spectra have been derived previously, the pseudo-continuum of 
the observed spectrum can be calculated through the analogous synthetic spectrum multiplied 
by a polynomial. In order to obtain a continuum fit over the whole range, the wavelength range 
is divided into four pieces (< 5600 A, 5400 - 6400 A, 6200 - 8000 A, > 7800 A) and a 5th 
order polynomial is applied to each of these ranges. After this piece-wise fitting, the pseudo¬ 
continuum of the entire wavelength range is fitted by a 11th order polynomial, then the original 
ID spectrum is normalized and delivered for the second stage. We made experiments of this 
normalization technique on the well calibrated SDSS/SEGUE (York et ah, 2000; Yanny et al., 
2009) survey spectra, and found that the difference between the determined parameters given 
by ULySS before and after normalization is negligible. However, when we applied this method 
to the LAMOST observations the differences are noticeable. This difference is mainly due to 
the unstable relative flux calibration of the spectra mentioned in Sec. 4.2.2; sometimes the 
shape of the flux-calibrated spectrum is quite strange, for example showing very steep regions 
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Fig. 5 An example of normalization through the CFI grid. Left: The pseudo¬ 
continuum derived from the corresponding parameter grid. Right: The normalized 
flux of the blue wavelength range (in blue) and the synthetic grid spectrum (in red), 
showing that the spectrum has been well fit by the CFI method. 


or with several big fake humps or other spurious features, which obviously cannot be fitted by 
a fixed order polynomial. On the other hand, this simple intervening normalization procedure 
can potentially alter of shape of the original physical lines. So besides the g-band SNR, the 
difference in the ULySS parameters between first and second stage results, to some extent, 
helps reflect the spectral data quality. If the two stages of ULySS measurements are similar to 
each other, ULySS determinations given by the normalized spectra are provided in the final 
publicly released AFGK stellar parameter catalog. Figure 6 illustrates one fitting example of 
ULySS for LAMOST observations. 

For observations with no prior knowledge, in order to prevent the ULySS algorithm from 
becoming trapped in a local minimum, we use a grid of initial guesses sampling all of the 
parameter space. For example, take the nodes of the guessing grid as T e ff =[4000, 7000, 15000] 
K, logg=[1.8, 3.8] dex, and [Fe/H] = [-1, 0] dex as adopted in the work of Wu et al. (2011a). The 
final solution (absolute minimum) is the best from those obtained with different guesses. For 
high quality AFGK type spectra, statistically by using a wide guessing grid, the global solution 
could be almost always be successfully found, with only several cases where the solution was a 
local minimum occurring out of ten thousand fitted stars. Also, for some spectra, the solution 
simply expressed by three physical parameters may not be unique, x 2 maps together with the 
convergence maps tool of the ULySS package (see (Wu et ah, 2011b)) were used to visualize 
the parameter space and identify the degeneracies and local minima. We can see that if taking 
the above guessing grid, each observation has to be fitted 12 times. For a survey project with 
massive data flow, it is quite time consuming, which is why we use CFI in front of ULySS to 
get a single initial guess. This effectively helps reduce the computational needs required, saving 
time and storage space for LASP by using fewer starting guesses in ULySS. 

Though CFI is just used to give a coarse estimation, for most cases its results used for 
the initial guesses in ULySS are acceptable (T e ^ is most reliable). In some cases, due to its 
low accuracy related either to the algorithm (those parameters are not physical isolated, or the 
spectral SNR if poor, for example) or from poor parameter coverage of the synthetic spectral 
grid, some of the determinations given by CFI are not credible. These biases are obvious for 
very hot and cold stars, and also for wide regions in the logg and [Fe/H] spaces. In practice, 
if any of the CFI parameters reside around the edges of its grid parameter coverage, or in an 
unphysical region, or out of the grid, ULySS does not take its results, but instead either uses a 
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Fig. 6 Fits of a LAMOST stellar spectrum with the ULySS ELODIE interpolator. The 
top panel represents the flux distribution - the observation is in black, and the model is 
in blue. The light blue is the multiplicative polynomial. Red lines represent the masked 
NaD telluric lines which were not well calibrated in ELODIE. The bottom panel 
shows the residuals between the observation and the best-fitted interpolated spec¬ 
trum (O-M). This spectrum shows star specid=20111108M31_012N38_B109032_v2.6.4 
with g-band SNR=17.51. The CFLstagel parameters are: [T e ff,logg, [Fe/H]] = 
[5734.07 K, 4.084 dex, -0.61 dex]; CFI_stage2 results: [5690.11 K, 4.01 dex, -0.74 dex], 
ULySS_stagel results: [5546.84 K, 4.31 dex, -0.40dex, -45.0 km/s], and ULySS_stage2 
results: [5605.01 K, 4.37 dex, -0.34 dex, -42.5 km/s]. 


moderate guess like [7000 K, 3.0 dex, 0.0 dex] or roughly estimates the starting value for T g q 
from the subclass estimated by the ID pipeline. 

Note that for the LAMOST DR1 data release, if the star is processed by LASP (i.e., it has 
T e ff,logg, and [FeH] values), the final released RV is determined by ULySS. If the spectrum 
was not processed by LASP, the RV is that measured by the ID pipeline, and for late M-type 
stars, the RV is from a separate module (inside the ID pipeline) which is specially tailored for 
cooler stars (Guo et ah, 2015). 

4-4-4 Why iterate once 

As known from the normalization described in section 4.4.3, the quality of the calculated pseudo¬ 
continuum depends on the stellar atmospheric parameters estimated by ULySS. There are large 
discrepancies between the synthetic and observed spectrum if the stellar atmospheric parameters 
from ULySS are unreliable. Therefore, the pseudo-continuum of observed spectrum calculated 
based on the synthetic spectrum is unacceptable, and the stellar atmospheric parameters ob¬ 
tained by ULySS using normalized spectra and parameters using un-normalized spectra will 
not agree. Likewise, if the two stellar atmospheric parameter estimations agree, it can be con¬ 
cluded that the stellar atmospheric parameters are reliable and the pseudo-continuum of the 
observed spectrum calculated based on the synthetic spectrum is acceptable. The reasonable 
degree of approximation can be measured statistically. In brief, if the two results approximately 
agree with each other, the stellar atmospheric parameters estimated using normalized spectra 
are adopted. 

4-4-5 Errors on the Stellar Atmospheric Parameters of LASP 

Since reliable random errors for each wavelength element are unavailable for LAMOST spectra, 
like in the work of Wu et al. (2011b), when using ULySS, we assumed a constant noise at all 
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wavelength points. We estimated an upper limit to the internal errors on the derived parameters 
by assuming % 2 = 1. We performed the fit with an arbitrary value of SNR and re-scaled the 
internal computed errors returned by ULySS by multiplying them by x 2 . The computed internal 
errors are small, likely for reasons related to the internal degeneracies between the parameters. 
We used statistics from comparison of the common targets observed by both LAMOST and 
SDSS/SEGUE (Wu et ah, 2014a, b) to estimate the external error. We then rescaled the errors 
using the ratio of the differences between LAMOST and SDSS parameters to the formal errors, 
assuming that the uncertainties of each series are equivalent. Currently the external errors 
determined in this way released in the parameters catalog are larger than the real measurement 
errors. 

Internal Comparison 

In total, about 18% of the targets in LAMOST have been observed more than one time. These 
multiple exposures can be used to estimate the internal errors on the parameters. We selected 
111,000 targets with two observations to estimate the internal errors. Figure 7 demonstrates 
the internal error distribution (i.e., the difference between parameters measured for the same 
star from different spectra) including a comparison of T e ff,log#, [Fe/H], and RV. 
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Fig. 7 Histograms comparing differences in measured stellar parameters for pairs of 
observations of the same star. The panels show the internal consistency of T e ff, log#, 
[Fe/H], and RV for 111,000 targets. 


Another way to estimate the internal error is comparing the parameters calculated through 
normalized spectra and un-normalized spectra. Figure 8 shows the distribution of difference 
between the parameters measured using normalized and un-normalized spectra. 
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Fig. 8 As in Figure 7, but now showing the difference between parameters measured 
using normalized and un-normalized spectra. 


Table 2 External Comparison with DR1 


T e ff error 

log# error 

[Fe/HJ error 

rv error 

High resolution 
SDSS DR9 
Gao et al. 
LSP3 

-47±95 K 
-91T111 K 
11±110 K 
7±133 K 

0.03T0.25 dex 
0.16T0.22 dex 
0.03T0.19 dex 
0.02T0.22 dex 

-0.02T0.1 dex 
0.14T0.05 dex 
0.01T0.11 dex 
0.06± 0.14 dex 

-7.2±6.6 km/s 
-3.78±4.91 km/s 
1.2± 3.1 km/s 


4-4-7 External errors 

We compared the common targets that have both LAMOST DR1 LASP measurements and high 
resolution spectral results (Bruntt et ah, 2012; Thygesen et ah, 2012; Molenda-Zakowicz et ah, 
2013; Dong et ah, 2014). For this comparison, T e ^ and [Fe/H] come from the high resolution 
spectra, while log# comes from Kepler (Borucki et ah, 2010) asteroseismic estimation. The 
comparison between LAMOST DR1 and high-resolution data is shown in Figure 9, along with 
statistical fits characterizing the quality of agreement. Figure 10 shows a similar comparison 
for the common targets between LAMOST (LASP results) and SDSS DR9. 

There are several independent works comparing parameters of DR1 with other reliable 
external databases. For example, Gao et al. compare DR1 parameters with the PASTEL catalog. 
Results from LSP3 (Xiang et al., 2015), an independent parameter pipeline developed at PKU, 
shows very good consistency with LAMOST/LASP measurements. All of these comparisons are 
summarized in Table 2. 

Zhang et al. (private communication) has also compared the parameters given by DR1 
with those from APOGEE, as shown in Figure 11. Besides these, Zhang et al. (2015) has 
compared 499 radial velocities (from DR1) with RVs of the same targets from MMT+Hectospec, 
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Fig. 9 Comparison of the three parameters measured by LASP on LAMOST DR1 
spectra and results from high resolution (HR) spectra (or Kepler asteroseismology, 
in the case of log#) of the same stars. u\ means the ULySS determination from the 
original spectra, and U 2 means the ULySS determination from the normalized spectra. 
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Fig. 10 Comparison of the three parameters between LASP and SDSS DR9 . The 
convention is same as described in Figure 9. 


which has RV accuracy of ^2.5 km s 1 (Drout et ah, 2009, 2012). The result is U r (DRl)- 
U r (MMT)=—3.76 zb 6.65 km s -1 (see Zhang et al. 2015). 

5 DATA PRODUCTS 

LAMOST DR1 consists of 2,204,860 spectra of stars, galaxies, quasars, and other unknown type. 
The footprint on the sky is shown in Figure 12, and the distributions of stellar parameters for 
all stars in RV vs. [Fe/H] and log# vs. T e ^ are shown in Figure 13. DR1 includes the following: 
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Fig. 11 Comparison of stellar parameters for stars in common between LAMOST 
DR1 and APOGEE. 


1. Spectra. - In general, there are 2,204,860 flux- and wavelength-calibrated, sky-subtracted 
spectra released publicly in the DR1. It includes 1,944,329 stars, 12,082 galaxies, 5,017 quasars, 
and 243,355 unknown objects. These spectra cover the range 3690-9100 A with a resolution of 
1800. 

2. Spectroscopic Parameters Catalogs. - In this data release, four spectroscopic parameters 
catalogs are included: the DR1 general catalog, the A type stars catalog, the late A, F, G and 
K type stellar parameter catalog, and the M star catalog. All of these include 30 basic pieces of 
information about each spectrum; e.g., RA, Dec, SNR, magnitudes, spectral type, and redshift, 
among others. In addition, the A type stars catalog also publishes four Balmer line indices and 
four spectral line widths at 20% below the local continuum, the late A, F, G and K type stars 
catalog provides effective temperature, surface gravity, and metallicity, and the M star catalog 
additionally contains a flag to indicate if a star shows magnetic activity. 

5.1 FITS File Description 

5.1.1 FITS File Designation 

The file names of the released FITS files are in the form ‘spec-MMMMM-YYYY_spXX- 
FFF.fits’, where ‘MMMMM’ represents the local modified Julian day (MJD), ‘YYYY’ is 
the plan identity string (PLANID), ‘XX’ indicates the spectrograph number (between 1 and 
16), and ‘FFF’ shows the fiber number (between 1 and 255). In addition, we provide a 
LAMOST designation for an object via the DESIG keyword that is formed as ‘LAMOST 
JHHMMSS.ss+DDMMSS.ss’, where ‘HHMMSS.ss’ is right ascension in units of HMS, and 
‘-hDDMMSS.ss’ is declination in units of DMS. 
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Fig. 12 The footprints of LAMOST DR1 in equatorial candidates. Pilot survey plates 
are shown as purple circles, and the first year data are green points. 



Fig. 13 The distribution of stars in LAMOST DR1 in the planes of RV vs. [Fe/H] 
(left) and log# vs. T e ff (right). 


5.1.2 Primary FITS Header 

We reorganize keywords of the primary FITS header in DR1 and divided them into eight groups 
including mandatory keywords, file information keywords, telescope parameter keywords, ob¬ 
servation parameter keywords, spectrograph parameters keywords, weather condition keywords, 
data reduction parameters keywords and spectra analysis results. See details in the appendix. 

5.1.3 Primary data array 

The primary data array has five rows and NAXIS1 (a keyword explained previously) columns. 
The five rows are Flux, Inverse Variance, Continuum-Subtracted Flux, Andmask and Ormask 
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of each. It is noted that the ‘inverse variance’ of the uncertainties (one over sigma-squared) 
can be used to estimate signal-to-noise ratio of each pixel ( flux * (inverse variance) 0 ' 5 ) , and 
‘Andmask’/‘Ormask’ are mask flags of six quality situations for each pixel, respectively. The 
details of these two masks are listed in Table.9 in the appendix. 

5.2 LAMOST Catalogs 

Four LAMOST DR1 catalogs are released and available at the website http://drl. lamost. 
org/, including the LAMOST general catalog, the A type stars catalog, the late A, F, G 
and K type stars catalog and the M star catalog. The detailed information is listed in tables 
(Table. 11,Table. 12,Table. 13) for each catalog. 

The uncertainties of measurements in DR1 are not included since the error estimation is 
unreliable and needs to be improved; however, the uncertainties will be released in DR2. 

5.2.1 LAMOST general catalog 

The LAMOST general catalog includes not only 717,660 spectra from the LAMOST pilot 
survey, comprised of 648,820 stellar spectra, 2723 galaxy spectra, 621 quasar spectra, and 65496 
unknown spectra, but also 1,487,200 spectra from the LAMOST general survey, which includes 
1,295,586 stellar spectra, 9,359 galaxy spectra, 4,396 quasar spectra, and 177,859 unknown 
spectra. In this catalog, there are 1,221,538 spectra with SNR in the g-band larger than 10, 
1,839,808 spectra with SNR in i-band larger than 10, and 1,198,410 spectra with SNR in both 
g- and i-bands larger than 10. 

All parameters of this catalog are listed in Table .10. For galaxies and quasars, ‘z’ provides 
redshifts for them, and their ‘z’ is set to -9999.00000000 if their redshift is unable to be estimated. 
For stars, ‘rv’ is their heliocentric radial velocity, and the Tv 5 is set to -9999.00000000 if their 
rv is unable to be estimated. For objects with unknown type, their ‘z’ and Tv’ are both set to 
-9999.00000000. In addition, for a galaxy or a quasar, the value of ‘rv’ is empty, and for a star, 
the value of ‘z’ is null. 

On the website http://drl.lamost.org, we provide two format with the same contents 
for LAMOST general catalogs; one catalog is in .fits format, and the other is .csv format. 

5.2.2 A type stars catalog 

There are 100,073 A-type stellar spectra in this catalog both from the pilot and general survey. 
Table. 11 shows all fields of this catalog. The ‘Class’ field presents the two-dimensional spec¬ 
tral classification results, and nearly all A type stars have luminosity class provided by the 
LAMOST analysis pipeline. Balmer line indices are included in ‘HalphaJndice’, ‘HbetaJndice’, 
‘HgamaJndice’, and ‘HdeltaJndice’ fields, and ‘Halpha_D0.2’, ‘Hbeta_D0.2’, ‘Hgama_D0.2’, 
and ’Hdelta_D0.2’ are the widths at 20% below the local continuum of these four Balmer 
lines. However, these eight fields are temporarily empty, and will be released in September 
of this year. The .csv format A type stars catalog can be download from the website (http: 
//drl. lamost. org). 

5.2.3 AFGK high quality stellar parameter catalog 

AFGK high quality stellar parameter catalog includes 61,686 A type stars, 548,214 F type stars, 
253,275 G type stars and 222,229 K type stars. These spectra are selected with SNR^ > 6 from 
dark-night observations and SNR^ > 20 from bright nights. The effective temperatures (T e ff), 
surface gravities (logg) and metallicities (Fe/H) are determined by LASP. The A type stars 
in this catalog are those from the A type star catalog with high SNR spectra. These data are 
available from the website http: //drl. lamost. org. 
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5.2.4 M star catalog 

The M star catalog contains 121,522 spectra, and their parameters are listed in Table.13. The 
basic parameters are provided as other catalogs, particularly the magnetic activity flags ‘magact’ 
are included in this catalog. The flag will be assigned to be 1 if stars have magnetic activity, 0 
if no magnetic activity and -9999 if the magnetic activity could not be determined due to low 
SNR. The specific catalog for M dwarfs is given in (Guo et ah, 2015). The .csv format M star 
catalog can be obtained from the website http: //drl. lamost. org. 

5.3 Caveats 

5.3.1 Radial velocity error 

The DR1 pipeline adopted local Beijing Time rather than UTC during the correction from 
geocentric coordinates to heliocentric coordinates. Thus an error in radical velocity is intro¬ 
duced, with mean value of 0.36 km s -1 . In addition, the RV errors of some A type stars reach 
80 km s -1 because their wavelength scale was incorrectly calibrated with air wavelengths instead 
of vacuum wavelengths. These are listed on the DR1 website. 

5.3.2 Spectral line errors coming from sub-exposure combination 

The LAMOST spectra released in DR1 are combined from several sub-exposures. If one of the 
three sub-exposures is problematic, spectral lines may be erroneously scaled (clipped) during 
the B-spline fitting by using the wrong data points. There are 66,838 spectra affected in DR1, 
which are listed on DR1 website. 

5.3.3 Offset of very bright targets 

Some coordinates of bright objects are artificially shifted to avoid saturated exposures. However, 
we list the original coordinates in the DR1 website. 

5.3.4 Mistaken spectra due to fiber problem 

Problems were found for tens of thousands of spectra in DR1 that are due to fibers or positioners. 
These were discovered after the data release, and include several thousands of spectra that have 
stellar parameters T e ff, log g , and [Fe/H]. A separate list of these problematic spectra is provided 
on the DR1 website along with their IDs, and will be excluded in DR2. 

6 DATA ACCESS AND USE 

The LAMOST DR1 can be accessed at LAMOST Data Release web portal, http://drl. 
lamost. org, which is developed under the collaboration with Chinese Virtual Observatory 
(China-VO). Web form based database search and data retrieve are available for everyone. 
Advanced features are provided for registered users, for example batch search, bulk retrieve 
and personal space. 

For registered users, lists of identifiers or coordinates supported by CDS SIMBAD 1 and 
VizieR 2 , or lists of LAMOST identifiers are acceptable to upload and use to search the database 
in batch mode. SQL-like queries are also supported. Virtual Observatory (VO) tools, such as 
Aladin 3 and SpecView 4 , and tools developed specifically for LAMOST are integrated into the 

1 http://simbad.u-strasbg.fr/simbad/ 

2 http://vizier.u-strasbg.fr/viz-bin/VizieR 

3 http://aladin.u-strasbg.fr/aladin.gml 

4 http://www.stsci.edu/institute/software_hardware/specview/ 
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web portal. Search results and spectra can be displayed and analyzed on the fly. Large results 
and datasets can be retrieved in several ways, including FTP, wget, Cloud storage, and VO 
data access interfaces. 

Based on CDS X-Match Service 5 , cross identification with dozens of well-known catalogs is 
also provided as an advanced feature of the data release portal. 

According to the LAMOST data policy, any publication making use of LAMOST data 
should cite the standard acknowledgement as following ’Guoshoujing Telescope (the Large Sky 
Area Multi-Object Fiber Spectroscopic Telescope LAMOST) is a National Major Scientific 
Project built by the Chinese Academy of Sciences. Funding for the project has been provided 
by the National Development and Reform Commission. LAMOST is operated and managed by 
the National Astronomical Observatories, Chinese Academy of Sciences.’ 
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APPENDIX : TABLE LIST 

There are nine tables for FITS files, including: Mandatory keywords(Table.l), File 
Information Keywords (Table. 2), Telescope Parameter Keywords(Table.3), Observation 
Parameter Key words (Table. 4), Spectrograph Parameters Keywords (Table. 5), Weather 
Condition Keywords (Table. 6), Data Reduction Parameters Keywords (Table. 7), Spectra 
Analysis Results Keyword (Table. 8) and the significance of six bits of ‘Andmask’ and 
‘Ormask’(Table.9). These values must be written in fixed format. 

There are another four tables, the gerneral catalog(Table.lO), A-type star catalog (Table. 11), 
Late A F G K catalog(Table.!2), and M-type star catalog(Table.!3). 


Table .1 Mandatory Keywords 


Keywords 

Value 

Comment 

SIMPLE 

T 

Primary Header created by MWRFITS vl.8 

BITPIX 

-32 


NAXIS 

2 


NAXIS1 

3909 


NAXIS2 

5 


EXTEND 

T 

Extensions may be present 


Notes: 

1. SIMPLE It is required to be the first keyword in the primary header of all FITS file. The value 
field shall contain a logical constant with the value T if the file conforms to this standard. This 
keyword is mandatory for the primary header and is not permitted in extension headers. A value of F 
signifies that the file does not conform to this standard. 

2. BITPIX The value field shall contain an integer, and it shall specify the number of bits that 
represent a data value. A value of -32 represents IEEE single precision floating point. 

3. NAXISThe value field shall contain a non-negative integer no greater than 999, representing the 
number of axes in the associated data array. A value of zero signifies that no data follow the header in 
the HDU. 

4. NAXIS1 and NAXIS2The value field of these two indexed keywords shall contain a 

non-negative integer, representing the number of elements along axis n of a data array. The NAXIS1 
keyword represents the number of wavelength array, i.e., the column number of the primary data 
array, and the NAXIS2 keyword indicates the row number of the primary data array. 

5. EXTENDThe value field shall contain a logical value indicating whether the FITS file is allowed 
to contain conforming extensions following the primary HDU. This keyword may only appear in the 
primary header and must not appear in an extension header. If the value field is T then there may be 
conforming extensions in the FITS file following the primary HDU. This keyword is only advisory, so 
its presence with a value T does not require that the FITS file contains extensions, nor does the 
absence of this keyword necessarily imply that the file does not contain extensions [6]. 
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Table .2 File Information Keywords 


Keywords 

Value 

Comment 

FILENAME 

AUTHOR 

l spec-55859-F5902_sp01-001.fits’ 

‘LAMOST Pipeline’ 

Who compiled the information 

N_EXTEN 

1 

The extension number 

EXTENO 

ORIGIN 

‘Flux Inverse Subcontinuum Andmask Ormask’ 
‘NAOC-LAMOST’ 

Organization 

DATE 

‘2013-07-17T10:23:12’ 

(UTC) 


Notes: 

1. FILENAME The value field shall contain a character string giving the name of this FITS file. 
Take the ‘spec-55859-F5902_sp01-001.fits’ as an example, ‘55859 is the local modified Julian day, 
‘F5902 is the plan ID, ‘spOl is the spectrograph ID, and ‘001 is the Fiber ID. 

2. AUTHOR This keyword contains a string constant LAMOST Pipline, which represents the 
author who produce this file. 

3. NJEXTEN The value field shall contain an integer giving the extension number of a FITS file. 

4. EXTENO This keyword contains a string constant Flux Inverse Subcontinuum Andmask Ormask 
explaining each row of the primary data array in a primary HDU. 

5. ORIGIN This ORIGIN keyword contains a string constant NAOC-LAMOST, which indicates the 
Organization responsible for this FITS file. NAOC represents the abbreviation of National 
Astronomical Observatories, Chinese Academy of Sciences. 

6. DATEThe value field shall contain a character string giving the UTC time when this FITS hie is 
created. 


Table .3 Telescope Parameter Keywords 


Keywords 

Value 

Comment 

TELESCOP 

‘LAMOST’ 

CuoShouJing Telescope 

LONGITUD 

117.580 


degree 

Longitude of site 

LATITUDE 

40.3900 


degree 

Latitude of site 

FOCUS 

19964 


mml Telescope focus 

CAMPRO 

‘NEWCAM’ 

Camera program name 

CAMVER 

‘v2.0 ’ 

Camera program version 


Notes: 

1. TELESCOP This keyword contain a string constant ‘LAMOST’ giving the name of our telescope. 

2. LONGITUD The keyword contain a floating-point constant, which provide the longitude of 
Xinglong station where LAMOST is mounted on. 

3. LATITUDE The keyword contain a floating-point constant, which provide the latitude of 
Xinglong station. 

4. FOCUS The keyword gives the telescope focus, and its unit is millimeter. 

5. CAMPRO The value held contain a string constant ‘NEWCAM’, which shows the name of 
camera. 

6. CAMVER The value held contain a character string ‘v2.0’, which gives the present camera 
program version. 
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Table .4 Observation Parameter Keywords 


Keywords 

Value 

DATE-OBS 

‘2013-06-03T16:47:08.38’ 

DATE-BEG 

‘2013-06-04T00:13:52.0’ 

DATE-END 

‘2013-06-04T01:21:17.0’ 

MJD 

56447 

MJDLIST 

‘56447-56447’ 

MJMLIST 

‘81283693-81283731’ 

PLANID 

‘HD174640N254456F01’ 

RA 

266.506521500 

DEC 

23.7254148000 

DESIG 

‘LAMOST J174601.56+234331.4’ 

FIBERID 

80 

CELLJD 

‘E0818 ’ 

X.VALUE 

51.4330635340 

Y.VALUE 

703.039784897 

OBJNAME 

‘230388178985208’ 

OBJTYPE 

‘Star’ 

OBJSOURC 

‘LEGUE.LCH’ 

FIBERTYP 

‘Obj’ 

MAGTYPE 

‘gri’ 

MAGI 

17.94 

MAG2 

17.47 

MAG3 

17.24 

MAG4 

99.00 

MAG5 

99.00 

MAG6 

99.00 

MAG7 

99.00 

OBS.TYPE 

‘OBJ’ 

OBSCOMM 

‘Science ’ 

RADECSYS 

‘FK5’ 

EQUINOX 

2000.00 


Comment 

The observation median UTC 
The observation start local time 
The observation end local time 
Local Modified Julian Day 
Local Modified Julian Day list 
Local Modified Julian Minute list 
Plan ID in use 

degree] Right ascension of object 
degree] Declination of object 
Designation of LAMOST target 
Fiber ID of Object 
Fiber Unit ID on the focal plane 

X coordinate of object on the focal plane 
Y coordinate of object on the focal plane 
Name of object 

Object type from input catalog 
Name of input catalog 
Fiber type of object 
Magnitude type of object 


magnitude 

magnitude 

magnitude 

magnitude 

magnitude 

magnitude 

magnitude 


Magi of object 
Mag2 of object 
Mag3 of object 
Mag4 of object 
Mag5 of object 
Mag6 of object 
Mag7 of object 

FLAT, 


The type of target (OBJ 
Science, Test 
Equatorial coordinate system 
Equinox in years 


ARC or BIAS) 


Notes: 

1. DATE-OBS The value field shall contain a character string, which gives the median moment UTC of multiple exposures. 

2. DATE-BEG The value field shall contain a character string giving the observation start Beijing Time. 

3. DATE-END The value field shall contain a character string, which provide the observation end Beijing Time. 

4. MJD The value field shall a non-negative integer giving the local modified Julian day. 

5. MJDLIST The value field shall contain a character string, which shows a list of local modified Julian day of n times 
exposures. 

6. MJMLIST The value field shall contain a character string, which shows a list of local modified Julian minute of n times 
exposures. 

7. PLANID The value field shall contain a character string providing the plan ID of the target. 

8. RA The value field shall contain a non-negative real floating-point number, which gives the right ascension of target. 

9. DEC The value field shall contain a non-negative real floating-point number, which gives the declination of target. 

10. DESIG The value field shall contain a character string, which indicates the name of LAMOST target. Like the name of 
SDSS target, numbers after the character J and before + represents RA in unit of HMS , and numbers after the character + are 
DEC in unit of DMS. 

11. FIBER.ID The value field shall contain a non-negative integer between 1 and 250, which shows the fiber ID and shall be 
used together with the spectrograph ID. 

12. CELL.ID The value field shall contain a character string, which gives the fiber unit ID on the focal plane. LAMOST focal 
plane is divided into four quadrant named EFGH respectively, the first character of this keyword represents the quadrant 
number, the first two numbers after the first character is the row number in this quadrant, and the next two numbers is the 
column numbers. 

13. X.VALUE and Y.VALUE Their value field shall contain two real floating-point numbers, which give X and Y 
coordinates of target on the focal plane. 

14. OBJNAME The value field shall contain character string, giving the name ID of object that determined by the RA, DEC 
and HTM method[7]. 

15. OBJTYPE The value field shall contain a character string giving the class of objects in input catalogs. 

16. OBJSOURC The value field shall contain a character string which shows the name of organization or person who submit 
input catalog. 

17. FIBERTYP The value field shall contain a character string, giving the type of fiber assigned to this target. This keyword 
has five values, i.e., Obj, Sky, F-std, Unused, PosErr and Dead. Obj means the fiber is assigned to a object, including star, 
galaxy and so on. Sky indicates the fiber is allocated to take skylight. F-std shows the fiber is used to take the light of a flux 
calibration standard star. Unused, PosErr and Dead mean the fiber is not used, goes to a wrong position, or does not work 
respectively. 

18. MAGTYPE The value field shall contain a character string, which shows the magnitude type of a target. 

19. MAGI, MAG2, MAG3, MAG4, MAG5, MAG6 and MAG7 The value field shall contain a real floating-point 
number between 0 and 100, giving the associated magnitudes of MAGTYPE keyword. For example, The MAGTYPE keyword is 
ugrizjh, the MAGI, MAG2, MAG3, MAG4, MAG5, MAG6 and MAG7 keywords provide the magnitudes of u, g, r, i, z, j and h 
filter respectively. 

20. OBS_TYPE The value field shall contain a character string giving the type of observation targets, which include object, 
flat, bias and arc lamp. 

21. OBSCOMMThe value field shall contain a character string constant representing the observation purposes , which 
includes observations used for science researches and kinds of tests. 

22. RADECSYSThe value field shall contain a character string giving the equatorial coordinate system based on the J2000 
position. 

23. EQUINOX The value field shall contain a real floating-point number giving the standard epoch used at present. 
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Table .5 Spectrograph Parameters Keywords 


Keywords 

Value < 

Comment 

SPID 

SPRA 

SPDEC 

SLIT_MOD 

1 ! 

266.640 
23.7144 
‘x2/3 ’ : 

Spectrograph ID 

degree] Average RA of this spectrograph 
degree] Average DEC of this spectrograph 
Slit mode, xl ,x2/3 or xl/2 


Notes: 

1. SPID The value field shall contain a non-integer numbers between 1 and 16, which provides the 
spectrograph ID. 

2. SPRA and SPDEC The value field of these two keywords shall contain two real floating-point 
numbers, which are the averages of RA and DEC of all objects in each spectrograph. 

3. SLIT_MOD The value field shall contain a character string giving the mode of slit, which 
includes xl, x2/3 and xl/2. At present, only mode x2/3 is available, which responds spectra 
resolution equals to 1800. 


Table .6 Weather Condition Keywords 


Keywords 

Value 

Comment 


TEMPCCDB 

-101.60 


C degree 

The temperature of blue CCD 

TEMPCCDR 

-103.90 


C degree 

The temperature of red CCD 

SEEING 

3.60 


arcsecond] Seeing during exposure 

MOONPHA 

24.51 


day] Moon phase for a 29.53 days period 

TEMP_AIR 

20.35 


C degree 

Temperature outside dome 

DEWPOINT 

3.11 


C degree 


DUST 

’ ’ 

Reservation 

HUMIDITY 

0.32 




WINDD 

87.00 


degree] Wind direction 

WINDS 

2.30 


m/s] Wind speed 


Notes: 

1. TEMPCCDB The value field shall contain a real floating-point number, which provides the 
temperature of blue CCD. The unit C degree represents centigrade degree. 

2. TEMPCCDR The value field shall contain a real floating-point number, which provides the 
temperature of red CCD. The unit C degree represents centigrade degree. 

3. SEEING The value field shall contain a real floating-point number giving seeing during exposure, 
which is calculated by manually measuring the full width at half maximum of guide star image. 

4. MOONPHA The value field shall contain a real floating-point number giving the moon phase. 

5. TEMP_AIR The value field shall contain a real floating-point number giving the temperature 
outside dome, which is measured by automatic weather instrument. The unit C degree represents 
centigrade degree. 

6. DEWPOINT The value field shall contain a real floating-point number giving the dew-point 
temperature, which is also measured by the automatic weather instrument. The unit C degree 
represents centigrade degree. 

7. DUST The value of this keyword is temporarily empty at present, because the dust measuring 
instrument is now in debugging, and we will write this parameters into fits header when problems are 
resolved. 

8. HUMIDITY The value field shall contain a real floating-point number between 0 and 1, which 
gives humidity in the air. 

9. WINDD The value field shall contain a real floating-point number which records the 
instantaneous wind direction when start exposure, and the direction of north is the 0 degree wind 
direction. 

10. WINDS The value field shall contain a real floating-point number which records the 
instantaneous wind speed when start exposure, and wind direction and speed are also measured also 
by the automatic weather instrument. 
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Table .7 Data Reduction Parameters Keywords(l) 

Keywords Value Comment 

LAMPLIST 4amphgcdne.dat’ Arc lamp emission line list 

SKYLIST ‘skylines.dat’ 


EXPID01 

‘01b-20130604001350-4-81283693’ 

EXPID02 

‘01b-20130604005112-5-81283731’ 

EXPID03 

‘01r-20130604001348-4-81283693’ 

EXPID04 

‘01r-20130604005109-5-81283731’ 

NEXP 

2 

NEXP.B 

2 

NEXP_R 

2 

EXPT_B 

3600.00 

EXPT_R 

3600.00 

EXPTIME 

3600.00 

BESTEXP 

81283693 

SCAMEAN 

2.11 

EXTRACT 

‘aperture’ 

SFLATTEN 

T 

PCASKYSB 

T 

NSKIES 

22 

SKYCHI2 

1.91072039307 

SCHI2MIN 

1.51493247984 

SCHI2MAX 

2.31929337901 


Sky emission line list 
ID string for exposure 1 
ID string for exposure 2 
ID string for exposure 3 
ID string for exposure 4 
Number of valid exposures 
Number of valid blue exposures 
Number of valid red exposures 
s] Blue exposure duration time 
s] Red exposure duration time 
s] Minimum of exposure time for all cameras 
MJM of the best exposure 
[ADU] Mean level of scatter light 
Extraction method 
Super flat has been applied 
PC A sky-subtraction has been applied 
Sky fiber number 
Meanchi 2 ofsky — subtraction 
Minimumchi 2 ofsky — subtraction 
M aximumchi 2 of sky — subtraction 


Notes: 

LAMPLIST The value field shall contain a character string giving the file name of arc lamp 
emission line list, which is used in the process of wavelength calibration. 

SKYLIST The value field shall contain a character string giving the file name of sky emission line 
list, which is used in the process of sky subtraction. 

EXPID01, EXPID02, EXPID03, EXPID04 The value field of these four keywords contain four 
character strings giving the ID strings of first, second, third and fourth exposures respectively, which 
are analogue to the file names of LAMOST raw data of four exposures. It should be noted that, if the 
exposure times of a target are n, there should be n associated keywords here, they are EXPID01, 
EXPID02, ....EXPIDn respectively. 

NEXP, NEXP_B, NEXP_R The value field of these three keywords shall contain three 
non-negative integers, which provide numbers of exposures, and numbers of valid blue and red 
exposures respectively. 

EXPT_B, EXPT_R The value fields of these two keywords shall contain two real floating-point 
numbers, which give exposure duration times of blue and red CCD. 

EXP TIME The value field shall contain a real floating-point, which gives the minimum of blue and 
red total exposures times. 

BESTEXP The value field shall contain a integer, which gives the MJM of a exposure with 
maximum signal and noise ratio in n time exposures. 

S CAME AN The value field shall contain a real floating-point giving the mean level of scatter light, 
which is the average flux of regions where there is no fiber and is at the left and right edge of a two 
dimension spectra image. 

EXTRACT The value field shall contain a character string, which indicates the method of spectrum 
extraction. In LAMOST spectra reduction pipeline, only the aperture method is applied to spectra 
extraction. 

SFLATTEN The value of this keyword shall be Boolean, which represents whether or not use the 
super flat. In LAMOST spectra reduction pipeline, super flat is used to make the fiber-to-fiber 
relative efficiency around 1, and it can be estimated through spline fitting the flat flux of all fibers in 
a spectrograph. 

PCASKYSB The value of this keyword shall be Boolean, which represents whether or not use the 
PCA method to subtract sky light. In LAMOST spectra reduction pipeline, the PCA method is used 
to subtract sky light at the wavelength range larger than 7200 angstrom. 

NSKIES The value field shall contain a integer, which shows the number of sky fiber in a 
spectrograph. 

SKYCHI2 The value field shall contain a real floating-point, which gives the mean chi square of 
sky-subtraction. In the process of LAMOST spectra reduction, super sky is obtained by spline fitting 
m sky spectra. And thus, the chi square between the super sky and each sky spectra in an exposure, 
and the average chi square of m sky spectra can also be able to obtain. Assuming n times exposures, 
there will be 2n average chi square because of n blue spectra and n red spectra, and this keyword will 
be evaluated by calculating the mean value of these 2n average chi squares. 

SCHI2MIN The value field shall contain a real floating-point, which gives the minimum chi square 
of sky-subtraction. As mentioned above, there will be 2n average chi squares assuming n time 
exposures, and this keyword will be the maximum of these chi squares. 

SCHI2MAX The value field shall contain a real floating-point, which gives the maximum chi square 
of sky-subtraction. As mentioned above, there will be 2n average chi squares assuming n time 
exposures, this keyword will be the minimum of these chi squares. 
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Table .7 Data Reduction Parameters Keywords(2) 


Keywords 

Value 

Comment 

NSTD 

2 

Number of (good) standard stars 

FSTAR 

‘149-182 ’ 

Fiber ID of flux standard stars 

FCBY 

‘catalog ’ 

Standard stars origin (auto, manual or catalog) 

HELIO 

T 

Heliocentric correction 

HELIO _RV 

-3.99846411542 

[km/s] Heliocentric correction 

VACUUM 

T 

Wavelengths are in vacuum 

NWORDER 

2 

Number of linear-loglO coefficients 

WFITTYPE 

‘LOG-LINEAR’ 

Linear-loglO dispersion 

COEFFO 

3.56820 

Central wavelength (loglO) of first pixel 

COEFF1 

0.000100000 

LoglO dispersion per pixel 

WATO-OOl 

‘system=linear’ 


WATU001 

wtype=linear label=Wavelength 



units=Angstroms 


CRVAL1 

3.56820 

Central wavelength (loglO) of first pixel 

CD1_1 

0.000100000 

LoglO dispersion per pixel 

CRPIX1 

1 

Starting pixel (1-indexed) 

CTYPE1 

‘LINEAR’ 


DC-FLAG 

1 

Log-linear flag 


Notes: 

NSTD The value field shall contain a non-negative integer, which shows the number of flux standard 
stars with good spectra quality. 

FSTAR The value field shall contain a character string giving the fiber identity numbers of flux 
standard stars, which are separated by the symbol 

FCBY The value field shall contain a character string giving the selection methods of flux standard 
stars, which include auto, manual and catalog. Auto represents the standard stars are selected by the 
LAMOST reduction pipeline, manual means they are picked out by experienced staffs, and catalog 
indicates the standard stars are provided by the input catalog. 

HELIO The value of this keyword shall be Boolean, which represents whether or not to perform the 
heliocentric correction. 

HELIO_RV The value field shall contain a real floating-point, which gives the radial velocity used to 
carry out the heliocentric correction. 

VACUUM The value of this keyword shall be Boolean, which represents whether or not the 
LAMOST spectra is converted to vacuum wavelength. 

NWORDER The value of this keyword shall contain a integer, which gives number of linear-loglO 
coefficients. 

WFITTYPE The value field shall contain a character string giving linear-loglO dispersion. 
COEFFO The value field shall contain a real floating-point number, which provides central 
wavelength (loglO) of first pixel. 

COEFF1 The value field shall contain a real floating-point number giving loglO dispersion per pixel. 
WATCLOOl The value field contain a character string. 

WAT1_001 The value field contain a character string. 

CRVAL1 The value field shall contain a real floating-point number, which gives the coordinate value 
of the reference pixel provided by the CRPIX1 keyword[8]. 

CD1_1 The value field shall contain a real floating-point giving the dispersion of per pixel. 

CRPIX1 The value of this keyword shall contain a integer, which sets the reference pixel location on 
pixel axis [8]. 

CTYPE1 The value field shall contain a character string, which will have the value LINEAR to 
define the wavelength axes to be linear [9]. 

DC-FLAG The value of this keyword shall be Boolean, a value of 0 defines a linear sampling of the 
dispersion and a value of 1 defines a logarithmic sampling of the dispersion[9]. 
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Table .8 Spectra Analysis Results Keyword 


Keywords 

Value 

Comment 

VERSPIPE 

V2.6.4 ’ 

Version of Pipeline 

CLASS 

‘STAR ’ 

Class of object 

SUBCLASS 

‘F5 ’ 

Subclass of object 

Z 

-0.000314280000000 

Redshift of object 

Z_ERR 

4 ’ 

Redshift error of object 

SN_U 

1.26 

SNR of u filter 

SN_G 

4.68 

SNR of g filter 

SN_R 

9.62 

SNR of r filter 

SNJ 

11.24 

SNR of i filter 

SN_Z 

7.44 

SNR of z filter 


Notes: 

VERSPIPE The value field shall contain a character string constant, which provides the version of 
LAMOST pipelines and are used to spectra processing and analysis. In this data release, the value of 
VERSPIPE is ‘v2.6.4. It should be noted that, ‘v2.6 is the version of spectra reduction pipeline, ‘v4 is 
the version of spectra analysis pipeline, and ‘v2.6.4 combines these two versions together. 

CLASS The value field shall contain a character string providing the classification result determined 
by the LAMOST spectra analysis pipeline, which includes ‘STAR, ‘GALAXY, ‘QSO or ‘Unknown. 
SUBCLASS The value field shall contain a character string, which gives a sub-classification results 
for stars. This keyword provides a more detailed spectra type for late A, F, G, K and M, and spectra 
and photometric type for A type stars. For galaxies, quasars or unknown type objects, this keyword is 
set to ‘Non now. 

Z The value field shall contain a real floating-point number providing redshift for a target, which is 
determined mainly by the LAMOST spectra analysis pipeline. For the case that redshift is unable to 
calculate by the pipeline, it will be manually determined through measuring the shifts of some 
spectral line centers. If the quality of a spectrum is poor, or it is classified as ‘Unknown, its redshift is 
artificially set to -9999.0. 

Z_ERR The value field shall contain a real floating-point number, which gives redshift error of a 
target. At present, the values of Z_ERR keywords are not yet published, and we will release redshift 
errors at the next data release in September of this year. 

SN_U, SN_G, SN_R, SN_I and SN_Z Keyword The value fields of these five keywords shall 
contain five real floating-point numbers, which give the signal and noise ratio (SNR) of u, g, r, i and z 
bands. Using the center wavelength and band width, we can obtain the wavelength range of each 
SDSS band, and then the SNR in each band is the median value at each pixel in this band. 


Table .9 The significance of six bits of ‘Andmask’ and ‘Ormask’ 


Bit 

Keyword 

Comments 

1 

BADCCD 

bad pixel on CCD 

2 

BADPROFILE 

bad profile in extraction 

3 

NOSKY no 

sky information at this wavelength 

4 

BRIGHTSKY 

sky level too high 

5 

BADCENTER 

fiber trace out of the CCD 

6 

NODATA 

no good data 


Notes: 

1. The ‘andmask’ information is a decimal integer determined by six-bit binary number, which 
represent six situations respectively as listed in this table. The associated bit of ‘andmask’ will be set 
to 1, if the case always appears in each exposure. Like ‘andmask’, ‘ormask’ information is also a 
decimal integer determined by six-bit binary number. The difference is that each bit of ‘ormask’ will 
be set to 1 if the related case happens in any exposure. 

2. For a spectrum, if you want to check in which case has happened during reduction process, you 
should convert the decimal ‘Andmask and ‘Ormaskto six-bit binary number, and a case must has 
happened in all exposures(at least one exposure) if associated bit is 1 in binary ‘Andmask(‘Ormask’). 
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Table .10 LAMOST general catalog 


Field (unit) 

Type 

Comment 

designation 

varchar 

Target Designation 

specObjld 

long integer 

Unique Spectra ID 

obsdate 

char 

Target Observation Date 

mjd 

char 

Local Modified Julian Day 

planid 

char 

Plan ID 

spid 

integer 

Spectrograph ID 

fiberid 

integer 

Fiber ID 

objra (degree) 

float 

Right Ascension 

objdec (degree) 

float 

Declination 

dataversion 

char 

Data Version 

snru 

float 

Signal Noise Ratio of u filter 

snrg 

float 

Signal Noise Ratio of g filter 

snrr 

float 

Signal Noise Ratio of r filter 

snri 

float 

Signal Noise Ratio of i filter 

snrz 

float 

Signal Noise Ratio of z filter 

objtype 

varchar 

Object Type 

class 

varchar 

Spectra Type 

subclass 

varchar 

Stellar Sub-Class 

magtype 

varchar 

Target Magnitude Type 

magi 

(mag) float 

Associated Magnitude 1 

mag2 

(mag) float 

Associated Magnitude 2 

mag3 

(mag) float 

Associated Magnitude 3 

mag4 

(mag) float 

Associated Magnitude 4 

mag5 

(mag) float 

Associated Magnitude 5 

mag6 

(mag) float 

Associated Magnitude 6 

mag7 

(mag) float 

Associated Magnitude 7 

tsource 

varchar 

Organization or person who submit input catalog 

tfrom 

varchar 

Input catalog submitted by an organization or a person 

fibertype 

varchar 

Fiber Type of target [Obj, Sky, F-std, Unused, PosErr, Dead] 

z 

float 

Redshift 

rv 

(km/s)float 

Heliocentric Radial Velocity 
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Table .11 A type stars catalog 


Field (unit) 

Type 

Comment 

designation 

varchar 

Target Designation 

specObjld 

long integer 

Unique Spectra ID 

obsdate 

char 

Target Observation Date 

mjd 

char 

Local Modified Julian Day 

planid 

char 

Plan ID 

spid 

integer 

Spectrograph ID 

fiberid 

integer 

Fiber ID 

objra (degree) 

float 

Right Ascension 

objdec (degree) 

float 

Declination 

dataversion 

char 

Data Version 

snru 

float 

Signal Noise Ratio of u filter 

snrg 

float 

Signal Noise Ratio of g filter 

snrr 

float 

Signal Noise Ratio of r filter 

snri 

float 

Signal Noise Ratio of i filter 

snrz 

float 

Signal Noise Ratio of z filter 

objtype 

varchar 

Object Type 

class 

varchar 

Stellar Class 

subclass 

varchar 

Stellar Sub-Class 

magtype 

varchar 

Target Magnitude Type 

magi (mag) 

float 

Associated Magnitude 1 

mag2 (mag) 

float 

Associated Magnitude 2 

mag3 (mag) 

float 

Associated Magnitude 3 

mag4 (mag) 

float 

Associated Magnitude 4 

mag5 (mag) 

float 

Associated Magnitude 5 

mag6 (mag) 

float 

Associated Magnitude 6 

mag7 (mag) 

float 

Associated Magnitude 7 

tsource 

varchar 

Organization or person who submit input catalog 

tfrom 

varchar 

Input catalog submitted by an organization or a person 

fibertype 

varchar 

Fiber Type of target [Obj, Sky, F-std, Unused, PosErr, Dead] 

rv (km/s) 

float 

Heliocentric Radial Velocity 

HalphaTndice (angstrom) 

float 

Line Indice of Halpha line 

HbetaJndice (angstrom) 

float 

Line Indice of Hbeta line 

HgamaTndice (angstrom) 

float 

Line Indice of Hgama line 

HdeltaJndice (angstrom) 

float 

Line Indice of Hdelta line 

Halpha_D0.2 (angstrom) 

float 

Width at 20% below the local continuum of Halpha line 

Hbeta_D0.2 (angstrom) 

float 

Width at 20% below the local continuum of Hbeta line 

Hgama_D0.2 (angstrom) 

float 

Width at 20% below the local continuum of Hgama line 

Hdelta_D0.2 (angstrom) 

float 

Width at 20% below the local continuum of Hdelta line 
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Table .12 A, F, G and K type stars catalog 


Field (unit) 

Type 

Comment 

designation 

varchar 

Target Designation 

specObjld 

long integer 

Unique Spectra ID 

obsdate 

char 

Target Observation Date 

mjd 

char 

Local Modified Julian Day 

planid 

char 

Plan ID 

spid 

integer 

Spectrograph ID 

fiberid 

integer 

Fiber ID 

objra (degree) 

float 

Right Ascension 

objdec (degree) 

float 

Declination 

dataversion 

char 

Data Version 

snru 

float 

Signal Noise Ratio of u filter 

snrg 

float 

Signal Noise Ratio of g filter 

snrr 

float 

Signal Noise Ratio of r filter 

snri 

float 

Signal Noise Ratio of i filter 

snrz 

float 

Signal Noise Ratio of z filter 

objtype 

varchar 

Object Type 

class 

varchar 

Stellar Class 

subclass 

varchar 

Stellar Sub-Class 

magtype 

varchar 

Target Magnitude Type 

magi (mag) 

float 

Associated Magnitude 1 

mag2 (mag) 

float 

Associated Magnitude 2 

mag3 (mag) 

float 

Associated Magnitude 3 

mag4 (mag) 

float 

Associated Magnitude 4 

mag5 (mag) 

float 

Associated Magnitude 5 

mag6 (mag) 

float 

Associated Magnitude 6 

mag7 (mag) 

float 

Associated Magnitude 7 

tsource 

varchar 

Organization or person who submit input catalog 

tfrom 

varchar 

Input catalog submitted by an organization or a person 

fibertype 

varchar 

Fiber Type of target [Obj, Sky, F-std, Unused, PosErr, Dead] 

rv (km/s) 

float 

Heliocentric Radial Velocity 

teff (K) 

float 

Effective Temperature 

log (dex) 

float 

Surface Gravity 

feh (dex) 

float 

Metallicity 



Data release one (DR1) of LAMOST general survey 
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Table .13 M-type star catalog 


Field (unit) 

Type 

Comment 

designation 

varchar 

Target Designation 

specObjld 

long integer 

Unique Spectra ID 

obsdate 

float 

Target Observation Date 

mjd 

char 

Local Modified Julian Day 

planid 

char 

Plan ID 

spid 

integer 

Spectrograph ID 

fiberid 

integer 

Fiber ID 

objra (degree) 

float 

Right Ascension 

objdec (degree) 

float 

Declination 

dataversion 

char 

Data Version 

snru 

float 

Signal Noise Ratio of u filter 

snrg 

float 

Signal Noise Ratio of g filter 

snrr 

float 

Signal Noise Ratio of r filter 

snri 

float 

Signal Noise Ratio of i filter 

snrz 

float 

Signal Noise Ratio of z filter 

objtype 

varchar 

Object Type 

class 

varchar 

Stellar Class 

subclass 

varchar 

Stellar Sub-Class 

magtype 

varchar 

Target Magnitude Type 

magi (mag) 

float 

Associated Magnitude 1 

mag2 (mag) 

float 

Associated Magnitude 2 

mag3 (mag) 

float 

Associated Magnitude 3 

mag4 (mag) 

float 

Associated Magnitude 4 

mag5 (mag) 

float 

Associated Magnitude 5 

mag6 (mag) 

float 

Associated Magnitude 6 

mag7 (mag) 

float 

Associated Magnitude 7 

tsource 

varchar 

Organization or person who submit input catalog 

tfrom 

varchar 

Input catalog submitted by an organization or a person 

fibertype 

varchar 

Fiber Type of target [Obj, Sky, F-std, Unused, PosErr, Dead] 

rv (km/s) 

float 

Heliocentric Radial Velocity 

magact 

integer 

Magnetic Activity 



